Strategies to mitigate shifts in red oak (Quercus sect. Lobatae) distribution under a changing climate.

Red oaks (Quercus sect. Lobatae) are a taxonomic group of hardwood trees, which occur in swamp forests, subtropical chaparral and savannahs from Columbia to Canada. They cover a wide range of ecological niches, and many species are thought to be able to cope with current trends in climate change. Genus Quercus encompasses ca. 500 species, of which ca. 80 make up sect. Lobatae. Species diversity is greatest within the southeastern USA and within the northern and eastern regions of Mexico. This review discusses the weak reproductive barriers between species of red oaks and the effects this has on speciation and niche range. Distribution and diversity have been shaped by drought adaptations common to the species of sect. Lobatae, which enable them to fill various xeric niches across the continent. Drought adaptive traits of this taxonomic group include deciduousness, deep tap roots, ring-porous xylem, regenerative stump sprouting, greater leaf thickness and smaller stomata. The complex interplay between these anatomical and morphological traits has given red oaks features of drought tolerance and avoidance. Here, we discuss physiological and genetic components of these adaptations to address how many species of sect. Lobatae reside within xeric sites and/or sustain normal metabolic function during drought. Although extensive drought adaptation appears to give sect. Lobatae a resilience to climate change, aging tree stands, oak life history traits and the current genetic structures place many red oak species at risk. Furthermore, oak decline, a complex interaction between abiotic and biotic agents, has severe effects on red oaks and is likely to accelerate species decline and fragmentation. We suggest that assisted migration can be used to avoid species fragmentation and increase climate change resilience of sect. Lobatae.

Genome-wide association study identifies SNP markers and putative candidate genes for terpene traits important for Leptocybe invasa resistance in Eucalyptus grandis.

Terpenes are an important group of plant specialized metabolites influencing, amongst other functions, defence mechanisms against pests. We used a genome-wide association study to identify single nucleotide polymorphism (SNP) markers and putative candidate genes for terpene traits. We tested 15,387 informative SNP markers derived from genotyping 416 Eucalyptus grandis individuals for association with 3 terpene traits, 1,8-cineole, γ-terpinene, and p-cymene. A multilocus mixed model analysis identified 21 SNP markers for 1,8-cineole on chromosomes 2, 4, 6, 7, 8, 9, 10, and 11, that individually explained 3.0%-8.4% and jointly 42.7% of the phenotypic variation. Association analysis of γ-terpinene found 32 significant SNP markers on chromosomes 1, 2, 4, 5, 6, 9, and 11, explaining 3.4-15.5% and jointly 54.5% of phenotypic variation. For p-cymene, 28 significant SNP markers were identified on chromosomes 1, 2, 3, 5, 6, 7, 10, and 11, explaining 3.4-16.1% of the phenotypic variation and jointly 46.9%. Our results show that variation underlying the 3 terpene traits is influenced by a few minor loci in combination with a few major effect loci, suggesting an oligogenic nature of the traits.

Characterization of terpene biosynthesis in Melaleuca quinquenervia and ecological consequences of terpene accumulation during myrtle rust infection.

Plants use a wide array of secondary metabolites including terpenes as defense against herbivore and pathogen attack, which can be constitutively expressed or induced. Here, we investigated aspects of the chemical and molecular basis of resistance against the exotic rust fungus Austropuccinia psidii in Melaleuca quinquenervia, with a focus on terpenes. Foliar terpenes of resistant and susceptible plants were quantified, and we assessed whether chemotypic variation contributed to resistance to infection by A. psidii. We found that chemotypes did not contribute to the resistance and susceptibility of M. quinquenervia. However, in one of the chemotypes (Chemotype 2), susceptible plants showed higher concentrations of several terpenes including α-pinene, limonene, 1,8-cineole, and viridiflorol compared with resistant plants. Transcriptome profiling of these plants showed that several TPS genes were strongly induced in response to infection by A. psidii. Functional characterization of these TPS showed them to be mono- and sesquiterpene synthases producing compounds including 1,8-cineole, ß-caryophyllene, viridiflorol and nerolidol. The expression of these TPS genes correlated with metabolite data in a susceptible plant. These results suggest the complexity of resistance mechanism regulated by M. quinquenervia and that modulation of terpenes may be one of the components that contribute to resistance against A. psidii.

A phylogenomic approach reveals a low somatic mutation rate in a long-lived plant.

Somatic mutations can have important effects on the life history, ecology, and evolution of plants, but the rate at which they accumulate is poorly understood and difficult to measure directly. Here, we develop a method to measure somatic mutations in individual plants and use it to estimate the somatic mutation rate in a large, long-lived, phenotypically mosaic Eucalyptus melliodora tree. Despite being 100 times larger than Arabidopsis, this tree has a per-generation mutation rate only ten times greater, which suggests that this species may have evolved mechanisms to reduce the mutation rate per unit of growth. This adds to a growing body of evidence that illuminates the correlated evolutionary shifts in mutation rate and life history in plants.

New approaches to tannin analysis of leaves can be used to explain in vitro biological activities associated with herbivore defence.

Although tannins have been an important focus of studies of plant-animal interactions, traditional tannin analyses cannot differentiate between the diversity of structures present in plants. This has limited our understanding of how different mixtures of these widespread secondary metabolites contribute to variation in biological activity. We used UPLC-MS/MS to determine the concentration and broad composition of tannins and polyphenols in 628 eucalypt (Eucalyptus, Corymbia and Angophora) samples, and related these to three in vitro functional measures believed to influence herbivore defence: protein precipitation capacity, oxidative activity at high pH and capacity to reduce in vitro nitrogen (N) digestibility. Protein precipitation capacity was most strongly correlated with concentrations of procyanidin subunits in proanthocyanidins (PAs), and late-eluting ellagitannins. Capacity to reduce in vitro N digestibility was affected most by the subunit composition and mean degree of polymerisation (mDP) of PAs. Finally, concentrations of ellagitannins and prodelphinidin subunits of PAs were the strongest determinants of oxidative activity. The results illustrate why measures of total tannins rarely correlate with animal feeding responses. However, they also confirm that the analytical techniques utilised here could allow researchers to understand how variation in tannins influence the ecology of individuals and populations of herbivores, and, ultimately, other ecosystem processes.

Overexpression of a developing xylem cDNA library in transgenic poplar generates high mutation rate specific to wood formation.

We investigated feasibility of the Full-length complementary DNA OvereXpression (FOX) system as a mutagenesis approach in poplar, using developing xylem tissue. The main goal was to assess the overall mutation rate and if the system will increase instances of mutants affected in traits linked to the xylem tissue. Indeed, we found a high mutation rate of 17.7%, whereas 80% of all mutants were significantly affected in cellulose, lignin and/or hemicellulose. Cell wall biosynthesis is a major process occurring during xylem development. Enrichment of mutants affected in cell wall composition suggests that the tissue source for the FOX library influenced the occurrence of mutants affected in a trait linked to this tissue. Additionally, we found that FLcDNAs from mutants affected in cell wall composition were homologous to genes known to be involved in cell wall biosynthesis and most recovered FLcDNAs corresponded to genes whose native expression was highest in xylem. We characterized in detail a mutant line with increased diameter. The phenotype was caused by a poplar homolog of LONELY GUY 1 (LOG1), which encodes an enzyme in cytokinin biosynthesis and significantly increased xylem proliferation. The causative role of LOG1 in the observed phenotype was further reaffirmed by elevated cytokinin concentration in the mutant and recapitulation overexpression experiment wherein multiple independent lines phenocopied the original FOX mutant. Our experiments show that the FOX approach can be efficiently used for gene discovery and molecular interrogation of traits specific to woody perennial growth and development.

Identifying genetic markers for a range of phylogenetic utility-From species to family level.

Resolving the phylogenetic relationships of closely related species using a small set of loci is challenging as sufficient information may not be captured from a limited sample of the genome. Relying on few loci can also be problematic when conflict between gene-trees arises from incomplete lineage sorting and/or ongoing hybridization, problems especially likely in recently diverged lineages. Here, we developed a method using limited genomic resources that allows identification of many low copy candidate loci from across the nuclear and chloroplast genomes, design probes for target capture and sequence the captured loci. To validate our method we present data from Eucalyptus and Melaleuca, two large and phylogenetically problematic genera within the Myrtaceae family. With one annotated genome, one transcriptome and two whole-genome shotgun sequences of one Eucalyptus and four Melaleuca species, respectively, we identified 212 loci representing 263 kbp for targeted sequence capture and sequencing. Of these, 209 were successfully tested from 47 samples across five related genera of Myrtaceae. The average percentage of reads mapped back to the reference was 57.6% with coverage of more than 20 reads per position across 83.5% of the data. The methods developed here should be applicable across a large range of taxa across all kingdoms. The core methods are very flexible, providing a platform for various genomic resource availabilities and are useful from shallow to deep phylogenies.

High marker density GWAS provides novel insights into the genomic architecture of terpene oil yield in Eucalyptus.

Terpenoid-based essential oils are economically important commodities, yet beyond their biosynthetic pathways, little is known about the genetic architecture of terpene oil yield from plants. Transport, storage, evaporative loss, transcriptional regulation and precursor competition may be important contributors to this complex trait. Here, we associate 2.39 million single nucleotide polymorphisms derived from shallow whole-genome sequencing of 468 Eucalyptus polybractea individuals with 12 traits related to the overall terpene yield, eight direct measures of terpene concentration and four biomass-related traits. Our results show that in addition to terpene biosynthesis, development of secretory cavities, where terpenes are both synthesized and stored, and transport of terpenes were important components of terpene yield. For sesquiterpene concentrations, the availability of precursors in the cytosol was important. Candidate terpene synthase genes for the production of 1,8-cineole and α-pinene, and ß-pinene (which comprised > 80% of the total terpenes) were functionally characterized as a 1,8-cineole synthase and a ß/α-pinene synthase. Our results provide novel insights into the genomic architecture of terpene yield and we provide candidate genes for breeding or engineering of crops for biofuels or the production of industrially valuable terpenes.

Accuracy of Genomic Prediction for Foliar Terpene Traits in Eucalyptus polybractea.

Unlike agricultural crops, most forest species have not had millennia of improvement through phenotypic selection, but can contribute energy and material resources and possibly help alleviate climate change. Yield gains similar to those achieved in agricultural crops over millennia could be made in forestry species with the use of genomic methods in a much shorter time frame. Here we compare various methods of genomic prediction for eight traits related to foliar terpene yield in Eucalyptus polybractea, a tree grown predominantly for the production of Eucalyptus oil. The genomic markers used in this study are derived from shallow whole genome sequencing of a population of 480 trees. We compare the traditional pedigree-based additive best linear unbiased predictors (ABLUP), genomic BLUP (GBLUP), BayesB genomic prediction model, and a form of GBLUP based on weighting markers according to their influence on traits (BLUP|GA). Predictive ability is assessed under varying marker densities of 10,000, 100,000 and 500,000 SNPs. Our results show that BayesB and BLUP|GA perform best across the eight traits. Predictive ability was higher for individual terpene traits, such as foliar α-pinene and 1,8-cineole concentration (0.59 and 0.73, respectively), than aggregate traits such as total foliar oil concentration (0.38). This is likely a function of the trait architecture and markers used. BLUP|GA was the best model for the two biomass related traits, height and 1 year change in height (0.25 and 0.19, respectively). Predictive ability increased with marker density for most traits, but with diminishing returns. The results of this study are a solid foundation for yield improvement of essential oil producing eucalypts. New markets such as biopolymers and terpene-derived biofuels could benefit from rapid yield increases in undomesticated oil-producing species.

Terpenes associated with resistance against the gall wasp, Leptocybe invasa, in Eucalyptus grandis.

Leptocybe invasa is an insect pest causing gall formation on oviposited shoot tips and leaves of Eucalyptus trees leading to leaf deformation, stunting, and death in severe cases. We previously observed different constitutive and induced terpenes, plant specialized metabolites that may act as attractants or repellents to insects, in a resistant and susceptible clone of Eucalyptus challenged with L. invasa. We tested the hypothesis that specific terpenes are associated with pest resistance in a Eucalyptus grandis half-sib population. Insect damage was scored over 2 infestation cycles, and leaves were harvested for near-infrared reflectance (NIR) and terpene measurements. We used Bayesian model averaging for terpene selection and obtained partial least squares NIR models to predict terpene content and L. invasa infestation damage. In our optimal model, 29% of the phenotypic variation could be explained by 7 terpenes, and the monoterpene combination, limonene, α-terpineol, and 1,8-cineole, could be predicted with an NIR prediction ability of  .67. Bayesian model averaging supported α-pinene, γ-terpinene, and iso-pinocarveol as important for predicting L. invasa infestation. Susceptibility was associated with increased γ-terpinene and α-pinene, which may act as a pest attractant, whereas reduced susceptibility was associated with iso-pinocarveol, which may act to recruit parasitoids or have direct toxic effects.

De Novo Transcriptome Study Identifies Candidate Genes Involved in Resistance to Austropuccinia psidii (Myrtle Rust) in Syzygium luehmannii (Riberry).

Austropuccinia psidii, causal agent of myrtle rust, was discovered in Australia in 2010 and has since become established on a wide range of species within the family Myrtaceae. Syzygium luehmannii, endemic to Australia, is an increasingly valuable berry crop. Plants were screened for responses to A. psidii inoculation, and specific resistance, in the form of localized necrosis, was determined in 29% of individuals. To understand the molecular basis underlying this response, mRNA was sequenced from leaf samples taken preinoculation, and at 24 and 48 h postinoculation, from four resistant and four susceptible plants. Analyses, based on de novo transcriptome assemblies for all plants, identified significant expression changes in resistant plants (438 transcripts) 48 h after pathogen exposure compared with susceptible plants (three transcripts). Most significantly up-regulated in resistant plants were gene homologs for transcription factors, receptor-like kinases, and enzymes involved in secondary metabolite pathways. A putative G-type lectin receptor-like kinase was exclusively expressed in resistant individuals and two transcripts incorporating toll/interleukin-1, nucleotide binding site, and leucine-rich repeat domains were up-regulated in resistant plants. The results of this study provide the first early gene expression profiles for a plant of the family Myrtaceae in response to the myrtle rust pathogen.

Transcriptome Profiling of Melaleuca quinquenervia Challenged by Myrtle Rust Reveals Differences in Defense Responses Among Resistant Individuals.

Plants have developed complex defense mechanisms to protect themselves against pathogens. A wide-host-range fungus, Austropuccinia psidii, which has caused severe damage to ecosystems and plantations worldwide, is a major threat to Australian ecosystems dominated by members of the family Myrtaceae. In particular, the east coast wetland foundation tree species Melaleuca quinquenervia, appears to be variably susceptible to this pathogen. Understanding the molecular basis of host resistance would enable better management of this rust disease. We identified resistant and susceptible individuals of M. quinquenervia and explored their differential gene expression in order to discover the molecular basis of resistance against A. psidii. Rust screening of germplasm showed a varying degree of response, with fully resistant to highly susceptible individuals. We used transcriptome profiling in samples collected before and at 5 days postinoculation (dpi). Differential gene expression analysis showed that numerous defense-related genes were induced in susceptible plants at 5 dpi. Mapping reads against the A. psidii genome showed that only susceptible plants contained fungal-derived transcripts. Resistant plants exhibited an overexpression of candidate A. psidii resistance-related genes such as receptor-like kinases, nucleotide-binding site leucine-rich repeat proteins, glutathione S-transferases, WRKY transcriptional regulators, and pathogenesis-related proteins. We identified large differences in the expression of defense-related genes among resistant individuals.

Temporal regulation of terpene synthase gene expression in Eucalyptus globulus leaves upon ozone and wounding stresses: relationships with stomatal ozone uptake and emission responses.

Ozone and wounding are key abiotic factors but, their interactive effects on temporal changes in terpene synthase gene expression and emission responses are poorly understood. Here, we applied combined acute ozone and wounding stresses to the constitutive isoprenoid-emitter Eucalyptus globulus and studied how isoprene, 1,8-cineole, and isoledene synthase genes were regulated, and how the gene expression was associated with temporal changes in photosynthetic characteristics, product emission rates, and stomatal ozone uptake through recovery phase. Photosynthetic characteristics and emission rate of isoprene, 1,8-cineole, and isoledene were synergistically altered, while three TPS gene expressions were antagonistically altered by combined stress applications. A time-delay analysis indicated that the best correspondences between gene expression and product emission rates were observed for 0 h time-shift for wounding and 0-2 h time-shifts for separate ozone, and combined ozone and wounding treatments. The best correspondence between ozone uptake and gene expression was observed for 0-4 h time-shifts for separate ozone and combined ozone and wounding treatments. Overall, this study demonstrated that expression profiles of isoprene, the monoterpene 1,8-cineole, and the sesquiterpene isoledene synthase genes differentially influenced their corresponding product emissions for separate and combined ozone and wounding treatments through recovery.

Complete mitochondrial genome of the green-lipped mussel, Perna canaliculus (Mollusca: Mytiloidea), from long nanopore sequencing reads.

We describe here the first complete genome assembly of the New Zealand green-lipped mussel, Perna canaliculus, mitochondrion. The assembly was performed de novo from a mix of long nanopore sequencing reads and short sequencing reads. The genome is 16,005 bp long. Comparison to other Mytiloidea mitochondrial genomes indicates important gene rearrangements in this family.

Four terpene synthases contribute to the generation of chemotypes in tea tree (Melaleuca alternifolia).

BACKGROUND: Terpene rich leaves are a characteristic of Myrtaceae. There is significant qualitative variation in the terpene profile of plants within a single species, which is observable as "chemotypes". Understanding the molecular basis of chemotypic variation will help explain how such variation is maintained in natural populations as well as allowing focussed breeding for those terpenes sought by industry. The leaves of the medicinal tea tree, Melaleuca alternifolia, are used to produce terpinen-4-ol rich tea tree oil, but there are six naturally occurring chemotypes; three cardinal chemotypes (dominated by terpinen-4-ol, terpinolene and 1,8-cineole, respectively) and three intermediates. It has been predicted that three distinct terpene synthases could be responsible for the maintenance of chemotypic variation in this species. RESULTS: We isolated and characterised the most abundant terpene synthases (TPSs) from the three cardinal chemotypes of M. alternifolia. Functional characterisation of these enzymes shows that they produce the dominant compounds in the foliar terpene profile of all six chemotypes. Using RNA-Seq, we investigated the expression of these and 24 additional putative terpene synthases in young leaves of all six chemotypes of M. alternifolia. CONCLUSIONS: Despite contributing to the variation patterns observed, variation in gene expression of the three TPS genes is not enough to explain all variation for the maintenance of chemotypes. Other candidate terpene synthases as well as other levels of regulation must also be involved. The results of this study provide novel insights into the complexity of terpene biosynthesis in natural populations of a non-model organism.

Genus-wide variation in foliar polyphenolics in eucalypts.

Many studies quantify total phenolics or total tannins, but understanding the ecological role of polyphenolic secondary metabolites requires at least an understanding of the diversity of phenolic groups present. We used UPLC-MS/MS to measure concentrations of different polyphenol groups - including the four most common tannin groups, the three most common flavonoid groups, and quinic acid derivatives - in foliage from 628 eucalypts from the genera Eucalyptus, Angophora and Corymbia. We also tested for phylogenetic signal in each of the phenolic groups. Many eucalypts contained high concentrations of polyphenols, particularly ellagitannins, which have been relatively poorly studied, but may possess strong oxidative activity. Because the biosynthetic pathways of many phenolic compounds share either precursors or enzymes, we found negative correlations between the concentrations of several of the constituents that we measured, including proanthocyanidins (PAs) and hydrolysable tannins (HTs), HTs and flavonol derivatives, and HTs and quinic acid derivatives. We observed moderate phylogenetic signal in all polyphenol constituents, apart from the concentration of the prodelphinidin subunit of PAs and the mean degree of polymerisation of PAs. These two traits, which have previously been shown to be important in determining plants' protein precipitation capacity, may have evolved under selection, perhaps in response to climate or herbivore pressure. Hence, the signature of evolutionary history appears to have been erased for these traits. This study is an important step in moving away from analysing "totals" to a better understanding of how phylogenetic effects influence phenolic composition, and how this in turn influences ecological processes.

Transcriptome analysis of terpene chemotypes of Melaleuca alternifolia across different tissues.

Plant chemotypes or chemical polymorphisms are defined by discrete variation in secondary metabolites within a species. This variation can have consequences for ecological interactions or the human use of plants. Understanding the molecular basis of chemotypic variation can help to explain how variation of plant secondary metabolites is controlled. We explored the transcriptomes of the 3 cardinal terpene chemotypes of Melaleuca alternifolia in young leaves, mature leaves, and stem and compared transcript abundance to variation in the constitutive profile of terpenes. Leaves from chemotype 1 plants (dominated by terpinen-4-ol) show a similar pattern of gene expression when compared to chemotype 5 plants (dominated by 1,8-cineole). Only terpene synthases in young leaves were differentially expressed between these chemotypes, supporting the idea that terpenes are mainly synthetized in young tissue. Chemotype 2 plants (dominated by terpinolene) show a greater degree of differential gene expression compared to the other chemotypes, which might be related to the isolation of plant populations that exhibit this chemotype and the possibility that the terpinolene synthase gene in M. alternifolia was derived by introgression from a closely related species, Melaleuca trichostachya. By using multivariate analyses, we were able to associate terpenes with candidate terpene synthases.

Identification of the Eucalyptus grandis chitinase gene family and expression characterization under different biotic stress challenges.

Eucalyptus grandis (W. Hill ex Maiden) is an Australian Myrtaceae tree grown for timber in many parts of the world and for which the annotated genome sequence is available. Known to be susceptible to a number of pests and diseases, E. grandis is a useful study organism for investigating defense responses in woody plants. Chitinases are widespread in plants and cleave glycosidic bonds of chitin, the major structural component of fungal cell walls and arthropod exoskeletons. They are encoded by an important class of genes known to be up-regulated in plants in response to pathogens. The current study identified 67 chitinase gene models from two families known as glycosyl hydrolase 18 and 19 (36 GH18 and 31 GH19) within the E. grandis genome assembly (v1.1), indicating a recent gene expansion. Sequences were aligned and analyzed as conforming to currently recognized plant chitinase classes (I-V). Unlike other woody species investigated to date, E. grandis has a single gene encoding a putative vacuolar targeted Class I chitinase. In response to Leptocybe invasa (Fisher & La Salle) (the eucalypt gall wasp) and Chrysoporthe austroafricana (Gryzenhout & M.J. Wingf. 2004) (causal agent of fungal stem canker), this Class IA chitinase is strongly up-regulated in both resistant and susceptible plants. Resistant plants, however, indicate greater constitutive expression and increased up-regulation than susceptible plants following fungal challenge. Up-regulation within fungal resistant clones was further confirmed with protein data. Clusters of putative chitinase genes, particularly on chromosomes 3 and 8, are significantly up-regulated in response to fungal challenge, while a cluster on chromosome 1 is significantly down-regulated in response to gall wasp. The results of this study show that the E. grandis genome has an expanded group of chitinase genes, compared with other plants. Despite this expansion, only a single Class I chitinase is present and this gene is highly up-regulated within diverse biotic stress conditions. Our research provides insight into a major class of defense genes within E. grandis and indicates the importance of the Class I chitinase.

Plant-Derived Terpenes: A Feedstock for Specialty Biofuels.

Research toward renewable and sustainable energy has identified specific terpenes capable of supplementing or replacing current petroleum-derived fuels. Despite being naturally produced and stored by many plants, there are few examples of commercial recovery of terpenes from plants because of low yields. Plant terpene biosynthesis is regulated at multiple levels, leading to wide variability in terpene content and chemistry. Advances in the plant molecular toolkit, including annotated genomes, high-throughput omics profiling, and genome editing, have begun to elucidate plant terpene metabolism, and such information is useful for bioengineering metabolic pathways for specific terpenes. We review here the status of terpenes as a specialty biofuel and discuss the potential of plants as a viable agronomic solution for future terpene-derived biofuels.

A curious case of resistance to a new encounter pathogen: myrtle rust in Australia.

Resistance genes (R genes) in plants mediate a highly specific response to microbial pathogens, often culminating in localized cell death. Such resistance is generally pathogen race specific and believed to be the result of evolutionary selection pressure. Where a host and pathogen do not share an evolutionary history, specific resistance is expected to be absent or rare. Puccinia psidii, the causal agent of myrtle rust, was recently introduced to Australia, a continent rich in myrtaceous taxa. Responses within species to this new pathogen range from full susceptibility to resistance. Using the myrtle rust case study, we examine models to account for the presence of resistance to new encounter pathogens, such as the retention of ancient R genes through prolonged 'trench warfare', pairing of resistance gene products and the guarding of host integrity.